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We put forward an approach to estimate the amount of bipartite spatial entanglement of down- 
converted photon states correlated in orbital angular momentum and the magnitude of the transverse 
(radial) wave vectors. Both degrees of freedom are properly considered in our framework which only 
requires azimuthal local linear optical transformations and mode selection analysis with two fiber 
detectors. The coincidence distributions predicted by our approach give an excellent fit to the 
distributions measured in a recent experiment aimed to show the very high-dimensional transverse 
entanglement of twin-photons from a down-conversion source. Our estimate for the Schmidt number 
is substantially lower but still confirms the presence of high-dimensional entanglement. 
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I. INTRODUCTION 

The phenomenon of quantum entanglement whereby 
distant systems can manifest perfectly random albeit per- 
fectly correlated behavior is now recognized as the es- 
sential ingredient to perform tasks which cannot be ac- 
complished with classically correlated systems [1 . The 
presence of entanglement has been traditionally revealed 
in the violation of Bell-type inequalities [2J. However, 
detecting such a violation does not provide in general a 
measure of the amount of entanglement. This is par- 
ticularly significant in systems correlated in multidimen- 
sional degrees of freedom [3J. Several techniques have 
been proposed to assess the presence of entanglement 
for different quantum scenarios. These include state to- 
mography (HG3E], that yields a complete reconstruction 
of a quantum state but requires many setting measure- 
ments, entanglement witnesses [7], which detect some 
entangled states with considerably less measurements, 
and the experimental determination of concurrence [HIS]- 
In Ref. [S], by using two copies of a down-converted 
two-photon state entangled in polarization and trans- 
verse momentum, the measurement of concurrence was 
achieved with a single, local measurement on one of the 
photons. 

Although bipartite entanglement is well understood, 
finding experimentally feasible procedures to quantify 
it for systems correlated in multidimensional degrees of 
freedom turns out to be quite challenging and relevant. 
Indeed, access to higher dimensional Hilbert spaces in 
which information can be encoded and manipulated has 
recently attracted great interest, with proof-of- principle 
demonstrations using quantum communication protocols 
in three-level systems (qutrits), such as entanglement 
concentration [10] , quantum bit commitment |llj and 
quantum coin-tossing |12j . Likewise, a complete char- 
acterization of states hyperentangled in polarization, or- 
bital angular momentum and frequency has been exper- 
imentally implemented |13j . 

The aim of this paper is to address the problem of 
how, by performing a certain set of local linear optical 
operations affecting one of two multidimensional spatial 



degrees of freedom in which two-photon states can be 
entangled, it is possible to obtain an explicit measure of 
the amount of bipartite transverse entanglement. Specif- 
ically, according to the interplay between both spatial 
degrees of freedom -orbital angular momentum (OAM) 
and the magnitude of the radial wave vectors- fundamen- 
tally different predictions are expected in the subsequent 
joint photodetection process, via mode-selection analysis 
with two fiber detectors preceded by azimuthal transfor- 
mations (acting only on the OAM). The application of 
our framework is compared with the results of a recent 
experiment |14j . 

The paper is organized as follows: Section II presents 
the general setting of the problem and includes the 
Schmidt decomposition technique for describing bipar- 
tite spatial entanglement. In Section III we reveal the 
generic features that arise in the photodetection coinci- 
dences according to the transverse structure of the two- 
photon wave function. Section IV illustrates the results 
found in Section III with an example of a realistic two- 
photon wave function that yields a full analytical solution 
to the problem. In Section V a closed-form expression 
for the Schmidt number is obtained in terms of easily 
accessible experimental parameters. We then proceed to 
exploit this Schmidt number in two interesting examples 
of optical transformations with azimuthal phase plates. 
These enable us to estimate the amount of spatial entan- 
glement of a two-photon source. Conclusions of the paper 
are drawn in Section VI. Details of our calculations, to- 
gether with some useful background material, have been 
included in two Appendices. 



II. MODAL SCHMIDT DECOMPOSITION FOR 
TWO-PHOTON STATES 

Two-photon pure quantum states are described in a 
Hilbert space by a continuous bilinear superposition of 
spatiotcmporal multimode states. Sources of such non- 
classical states of light are mostly realised in the process 
of spontaneous parametric down-conversion |15j . where 
an intense quasimonochromatic laser pump illuminates 
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a crystal endowed with a quadratic nonlinearity produc- 
ing pairs of photons (idler and signal). Conservation of 
energy and momentum impose that the state be spec- 
trally and spatially correlated. Here we explore entan- 
glement involving spatial degrees of freedom that de- 
pend on the transverse structure of these states. For 
simplicity, we assume that the down-converted photons 
are linearly polarized, monochromatic and frequency de- 
generated. The two-photon state can then be written as 
\ip) = J dq i dq s <J>(q i ,q s )a t (q i )a t (q ;s )|vac), where q i s are 
the transverse components of the idler and signal wave 
vectors. Under conditions of paraxial and nearly collincar 
propagation of the pump, idler and signal photons, the 
two-photon amplitude <& is given by [16] 4>(qj,q s ) — 
E((\i + q s )G(qi — q s ). The function E represents the 
transverse profile of the pump, whereas G originates from 
the phase-matching condition in the longitudinal direc- 
tion and depends on the specific orientation and cut of 
the nonlinear crystal. Since the arguments of E and G 
enforce correlations in different manifolds of the idler and 
signal wave vector space, it is the global structure of <& 
the one dictating the entanglement degree of \tp). 

In order to extract the amount of entanglement con- 
tained in one may resort to the Schmidt decomposi- 
tion that provides the spatial information modes of the 
two- photon pair. Suppose that the transverse spatial fre- 
quency field of the pump beam has a Gaussian profile of 
the form E(q.i + q s ) cx e - ™o' qi+qs '"/ 4 , where Wq is the 
width (at the beam waist). The chosen pump profile 
peaks when its argument vanishes. This imposes that 
the idler and signal transverse wave vectors should be 
mostly anticorrelated (q^ ~ — cb)- Remarkably, it was 
shown by Law and Eberly |17j that with such a Gaussian 
profile E the normalised two-photon amplitude can be 
expressed as 

oo oo 
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where U£,„(qi) = e I ^^ ! „(g i )/v / 27r and u_^ n (q s ) = 
e - ll <t's V _ l in (g s )/v2~7r are the normalised polar Schmidt 
mode functions for the idler and signal photons with 
topological charge I and radial index n corresponding 
to eigenvalues \i n (they satisfy 1 > \i n > and 
J2e n ^<"n = 1) of the reduced density matrices for each 
photon. Knowledge of Xg n yields a direct measure of the 
degree of transverse entanglement given by the Schmidt 
number [18] K, — (J2e n ^n) _1 > which is the reciprocal 
of the purity of the idler and signal density matrices, 
it is invariant under free propagation and yields an av- 
erage of the number of relevant spatial modes involved 
in the decomposition. The larger the value of /C, the 
higher the transverse entanglement. For instance, prod- 
uct states correspond to K — 1 (there is only one nonvan- 
ishing eigenvalue equal to 1), whereas states with JC > 1 
are entangled. A distinguishing feature of decomposi- 
tion is that it represents $ in terms of a perfectly 
correlated discrete basis of paraxial eigenstates of the 




FIG. 1: (color online). Two-photon coincidence detection 
configuration. Down-converted idler and signal beams from a 
nonlinear crystal traverse two optical systems (h^' s ') which 
include azimuthal phase plates and are coupled into single- 
mode fibers (SMF) gated by a coincidence circuit (CC). 

OAM operator along the direction of light propagation 
(with corresponding eigenvalues £h) [19] [20] , rather than 
in a continuous plane- wave modal expansion. The precise 
form of the radial idler and signal eigenmodes ve t n(<li) 
and v-£^ n (q s ) depends on the specific phase-matching 
function G. In particular, when G is approximated by 
a constant there are striking consequences: Eq. (JT|) be- 
comes a (non-normalisable) superposition in which all 
OAM eigenstates have an equal weight and the radial de- 
pendence becomes just a global factor. It is important to 
emphasize that to attain decomposition the appropri- 
ate choice of the widths, Wi and w s , for the idler and sig- 
nal radial eigenmodes, has to be made. If the two-photon 
amplitude is of the form $(qj,q s ) = _B(qj+q s )G(q i — q s ), 
then Wi = w s = ws, where w$ is the so-called Schmidt 
width. For widths different from an additional sum- 
mation over the radial indexes occurs, and the perfect 
correlation between the idler and signal radial modes is 
absent (see Appendix A). Moreover, the fact that the 
idler and signal mode functions are anticorrelated with 
respect to their topological charge numbers is a conse- 
quence of a more general process: the conservation of 
OAM, which is transferred from the pump photon (car- 
rying zero OAM for a Gaussian mode) to the down- 
converted photon pair 21J. 

III. AZIMUTHAL TRANSFORMATIONS ON 
THE TWO-PHOTON STATE 

The usefulness of the Schmidt decomposition Q be- 
comes apparent when analyzing the propagation of two- 
photon states through optical systems. Each of the in- 
tervening modes evolves and transforms independently of 
the others. It is also clear that the correlation proper- 
ties displayed by the two-photon amplitude are pre- 
served in the position representation. Therefore, sup- 
pose that the idler and signal photon beams, described 
now by the paraxial two-photon wave function in the 
transverse position representation E3] ip( r i, r s) — 

(ri,r s |V>) = Y,e,n(~ i y^^ u ^A r i) u -<i,n( r s), are each 
transmitted through different linear optical systems that 
include diffractive (or refractive) azimuthal phase plates 
(see Fig. [T|. The role of the plates in each path is 
to imprint an azimuth-dependent phase factor on the 
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incoming Schmidt modes. Their (separate) action on 
the two-photon wave function can be represented by 
the unitary and radially symmetric impulse response 
functions h^' s '((f),(j/). These functions locally trans- 
form the spiral harmonic mode content of each pho- 
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J o 2 "/ o 2 "d ( /.d0'e^e- 4£ >'/i(^ s )(0>')/27r, so that the re- 
sulting output two-photon wave function is 
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where the initial perfect correlation in OAM is lost, re- 
maining only that corresponding to the radial modes. 

Upon traversing their respective linear optical systems 
the idler and signal photons, having OAM £; — £a, £s — 
£b and radial indexes ua,b, are detected in coincidence. 
By placing photodetectors at the output ports of suitable 
arrays of deterministic mode sorter interferometers |24l 
[25] . or computer-generated holograms [TUj, it is possible 
to distinguish modes bearing different OAM. In practice, 
the most straightforward procedure involves projecting 
into single-mode fibers, where the propagated mode has a 
fundamental Gaussian profile (£a = £b = = jib = 0). 

The probability that the idler and signal photons will 
be projected into modes ui A nA and ug B nB is found to 
beP^? = \jdv l dv s ul An ^v i )u* lBnB {v s )^ out (v l ,v s )\\ 
with an additional incoherent (and weighted) sum over 
the radial indexes Ua,b when taking into account multi- 
mode detection. This photodection probability is ex- 
pressed in terms of the spatial overlap of the Schmidt and 
the fiber modes at the planes where the phase plates are 
located. Notice that their corresponding widths, u>s and 
wg, are not necessarily equal, and thus the orthogonal- 
ity between these modes when their radial indexes differ 
does not hold in general. To evaluate V™ A '™ B , we de- 
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xv e~n ' ( r i) v -en( r s)> wnere the two different widths of 
the radial modes have been specified for clarity. The co- 
incidence probability can then be written as 
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All the radial dependence is contained in the functions 



1Zf' B that modulate the angular impulse response func- 
tions h( l ' s K We emphasize that although the radial part 
of the Schmidt modes does not experience any signif- 
icant transformation when the idler and signal beams 
traverse their respective linear optical systems, its proper 
inclusion in the detection process is essential. This is re- 
flected in the structure of Eq. ([3| with the presence of 
^-dependent radial functions 1Z £ ' , and is a consequence 
of a non-constant phase-matching function G. Had G 
been a constant then the input two-photon wave func- 
tion would have been represented by a common radial 



function times a (non-normalisable) maximally entangled 
superposition of spiral harmonic modes. In this limit- 
ing case one derives a fundamentally different prediction 
for the coincidences: Ve A: e B 



where now all the radial dependence of the detected 
modes appears only as a global function. 



IV. AN EXACTLY SOLUBLE MODEL FOR THE 
TWO-PHOTON AMPLITUDE 

To illustrate the previous results, suppose that the 
phase-matching function G is also Gaussian (see Ap- 
pendix C, where a justification of this model is provided), 
an approximation implying that photons are generated 
near the phase-matching region of wave vectors where 
the down-conversion process occurs most effectively. The 
normalised two-photon amplitude reads 



(q . i;qs ) = ^ e -^|q.+q s | 2 /4 e - h 2 |q,-q 



a | 2 /4 



(4) 



where b < wo plays the role of an effective width of 
the phase-matching function, and depends on the non- 
linear crystal thickness, although no specific relation is 
assumed here. At variance with other models, where 
G is often approximated by a constant (b — > 0), the 
two-photon amplitude ^ captures the relevant features 
of the transverse wave-vector correlation between the 
idler and signal photons, and provides an analytically 
amenable model that yields explicit formulas for Eq. Q . 
We show in Appendix A that for the two-photon ampli- 
tude Q , the Schmidt mode functions belong to the well- 
known Laguerre-Gaussian basis [17 , or, more generally, 
to a continuum family of spatial modes generated via 
metaplectic mappings (rotations on the orbital Poincare 
sphere) from the Laguerre-Gaussian modes [20j [26j EI] . 
Furthermore, we also find that the Schmidt eigenvalues 
are \ ln = (1 - f)*?W+*>, with £ = (w - b)/(w + b). 
This allows us to express the Schmidt number in the fol- 
lowing closed form JC = [(1 + £ 2 )/(l — £ 2 )] 2 - Large values 
of /C occur when 6^0. The minimum K, = 1 is attained 
when b — > wq (the two-photon amplitude Q becomes a 
separable function in the idler and signal wave vectors). 
The Schmidt width for any of the above families of eigen- 
modes is always the same: ws = \/2wob. We stress that 
the Schmidt width does not represent the actual cross- 
section widths Wi.s of the idler and signal beams. These 
widths, that can be measured experimentally, can also 
be obtained by resorting to the partially reduced density 
matrix of the idler and signal photons. For amplitude Q 
one finds Wi lS = ^/2wq + b 2 , which is consistent with the 
widths employed in [28 in the thin crystal approximation 
(6-0). 

Let us focus on the most usual encountered situation 
where the measured fiber modes are the fundamental 
Gaussian modes (of width wg at the phase plates). Since 
the involved eigenfunctions in the Schmidt decomposi- 
tion have cylindrical symmetry, we use the Laguerre- 
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Gaussian modes as the convenient computational basis to 
derive the radial functions IZi = ji^ A ' B ~ nA - B ~ Q ( see Ap- 
pendix B). Remarkably, it turns out that the functions 
IZi can be cast in terms of a sole parameter s 



r*(i + f) \i\ 



\t\;s 2 ) s W 



(5) 



where F(a, b; c; d) denotes the hypergeometric function, 
and 

2£ 



i+e + a-e)(ws/w G r 

w 2 ~b 2 
w 2 + b 2 + (2w b/w G ) 2 ' 



(0) 



The parameter s corrects £ by taking into account the 
spatial overlapping of the Schmidt and the fiber modes. 
Indeed, if and only if ws = wg, does s — £ hold. When 
b — > then s — > 1, which corresponds to a constant 
phase-matching function. The functions IZi(s) increase 
monotonically from 72.^(0) = to TZe(l) = 1 for £ ^ 
(TZo(s) = 1). The fact that all IZe(s) depend on a single 
parameter s will be exploited in the next Section to show 
how one can estimate the Schmidt number in a particular 
experimental scenario. 



V. EXPERIMENTAL SCHMIDT NUMBER 

Let us now examine the problem of measuring the 
amount of transverse entanglement of two-photon sources 
which, in our case, is characterised by the Schmidt num- 
ber /C. In principle, a complete quantum tomography 
of the two-photon state could yield the desired amount 
of entanglement but generally this would require a 
very large number of measurements, each for every possi- 
ble pair of spatial modes, and this is technically demand- 
ing. Here we propose an alternative approach. In our 
simple model for the two-photon amplitude Q we have 
found that the radial part of the coincidence probabili- 
ties ([3]) depends on a single parameter s. This parame- 
ter s involves the characteristic widths appearing in the 
two-photon amplitude Q , namely the pump beam width 
Wq and the phase-matching width b, together with the 
fiber mode width wg (at the phase plate locations). The 
phase-matching width b could be measured by scanning 
in the plane of detection, but in our case it is not neces- 
sary. If, instead, one rotates the azimuthal phase plates 
(maintaining the detectors fixed), then, the recording of 
the coincidence distributions allows one to extract the 
value of s as a fitting parameter for V s = 
via Eqs. ^ and (J5j. Therefore, if s is conceived as a 
parameter to be directly measured (rather than b), the 
amount of transverse entanglement can now be written 
in the form 
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FIG. 2: (Left column) Normalised coincidence distributions 
as a function of the relative orientation between the idler and 
signal angular diaphragms. Solid, dashed and dotted curves 
are calculated for 8 = 1, 0.7, and 0.4, respectively, (a) f3 — n/2 
and (c) (3 — ir, both with n = 0.5. (Right column) The 
corresponding configurations of the angular diaphragms. 



where /1 = Wq/wg- This experimental Schmidt number 
depends on quantities that are easily accessible: the fit- 
ting parameter s and the widths u>o and wg (the presence 
of the width ratio /i should be interpreted as a correcting 
geometrical factor). It increases from K, = 1, when s = 0, 
to infinity as s — ► 1. 

By properly engineering the impulse response func- 
tions h^-^ it is possible to enhance the sensitivity of V s 
with s, thus improving the accuracy of the estimated IC. 
To this end, we consider two simple types of transpar- 
ent azimuthal phase plates, and examine the dependence 
of Eq. ([3| when the idler and signal phase plates are 
mutually rotated a relative angle a = ai — a s . Their 
dispersionless impulse response functions are of the form 
h^(<t),<j)') = e ie ^^5{ct)-(t)'). Owing to the initial per- 
fect anticorrelation in OAM of the down-converted pho- 
ton pairs, one expects that only when both photons are 
subjected to complementary azimuthal transformations 
the perfect anticorrelation is preserved and the coinci- 
dences are maximal. As soon as the phase plates are ro- 
tated in such a way that they are not longer oriented in 
a complementary arrangement, the photon coincidences 
decrease. Notice that the axes with respect of which the 
angles and a s are taken do not coincide. For symmet- 
ric phase plates these reference axes are inverted 180°. In 
what follows, we assume that each plate is characterised 
by a noninteger parameter rj — (n — l)d/\, where d is 
the relative step height introduced by the plates, no their 
refractive index, and A the wavelength of the idler and 
signal beams. 

The first phase plate type we consider is an an- 
gular diaphragm [see inset in Fig. [2ja)] . It con- 
sists of a thin uniform dielectric circular slab with 
a "cake-slice" indentation that subtends an angle (5, 
with a nonzero 6i(4>) = —6 s {(j)) = 2nr] (mod 2tt) 
only if ai tS + (3 < cj) < ai tS + 2ir. Similar angu- 
lar diaphragms have been employed in proof-of-principlc 
demonstrations of the uncertainty relation for angu- 
lar position and OAM 29J. The probability (|3| can 

be cast as T { s AD \a) = [Ko(s)[(0 - ir) 2 + ir 2 cot 2 (7rn)] 
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FIG. 3: Normalised coincidence distributions as a function 
of the relative orientation between the idler and signal spiral 
phase plates. Solid, dashed and dotted curves are calculated 
for s = 1,0.7, and 0.4, respectively, (a) rj = 0.5; and (b) 
rj = 4.5. (c) Configuration of the spiral phase plates. 



+8^fci^(s)cos(£a)sin 2 (£ ) 5/2)/^ 2 ] . The main fea- 
tures of the normalized coincidence vi AD \a) /vi AD \o) 
are: (i) it does not depend on the integer part of rj 
(it suffices to consider < 77 < 1), the visibility be- 
ing maximal when 77 = 1/2; (ii) the coincidence dis- 
tributions are identical whether the aperture angle is 
[3 or 2ir — [3 and symmetrical around a — tt; (iii) for 
P =/= 7r and rj fixed, the visibility diminishes as s decreases; 
and (iv) the maximum visibility always occurs in the 
limit s — > 1 (constant phase-matching function and thus 

very high Schmidt number) where one has V { £?(a) = 
tt 2 [7r[cot 2 (7r?7) - 1] + |27r - a- (3\ + \a- f3\] 2 . Figures 
[2ja) and (c) depict the characteristic profiles of the (nor- 
malised) vi AD ^ when the aperture angles of the angular 
diaphragms are (3 = n/2 and /3 = tt [Figs. [2jb) and 
(d) show, respectively, the configurations of the angu- 
lar diaphragms that yield the coincidences (a) and (c)]. 
Comparing Figs. |2ja) and (c) one sees that the visibility 
of the former exhibits a much stronger variation with s 
than the latter. This suggests that the phase plate con- 
figuration of Fig. |2jb) is preferred to that of Fig. [2^d) 
for achieving a more accurate estimation of the Schmidt 
number via Eq. 0. Finally, we should add that if rj is an 

integer number then the profiles of vi AD ^ are constant 
(independent of a). 

The second type of azimuthal (refractive) component 
is a spiral phase plate [see inset in Fig. [3]ja)] [T9j 150] . 
Its impulse response function includes the phase de- 
pendence 8i{4>) = —9 S {4>) = rj(f>, for cvi. s < (j) < 
ai tS + 2tt. In this case, the joint probability (J3J) can be 



written as V. 



(SPP) 



(a) = \J2 e K e (s)e^(t + rj)- 



At 



variance with vi AD \ the coincidence distributions ex- 
hibit the development of interference ripples as s de- 
creases and rj increases [see Figs. |3ja) and (b)]. When 

s — > 1 a parabolic profile is obtained (a) — 

tt 2 [tt 2 cot 2 (7r?7) + (a — n) 2 ] I sin 2 (7r?7). In this limit the 
coincidences become independent of the integer part of 
rj, and, as with angular diaphragms, the maximum coin- 
cidence visibility is attained (when rj = 1/2). 
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FIG. 4: Coincidence distribution dependence on the relative 
orientation between idler and signal spiral phase plates (SPP 
with rj — 3.5). Circles are experimental values from Ref. |14j . 
solid and dashed curves correspond to "pi SPP ' fitted with s = 
0.66, and Pg^[ P \ respectively. The slight asymmetry in the 
experimental values was probably caused by the presence of 
a small anomaly in the center of the spiral phase plates [14] . 
(Inset) Schmidt number A3 as a function of the width ratio 
ji — wo/wa for s = 0.66. 



We have considered other types of azimuthal phase 
plates and found the same parabolic dependence of the 
coincidences as s — > 1 . This is consistent with the above- 
mentioned prediction that for very high transverse en- 
tanglement the joint probability becomes independent of 
the radial structure of the Schmidt modes. On the other 
hand, the absence of a vanishing minimum in the co- 
incidence distributions is a signature of a finite amount 
of transverse entanglement. This effect gives rise to a 
contribution in the photocounts which always exists re- 
gardless of the detector efficiencies [see Fig. IJa) and 
Figs. [3ja),(b)]. It can however be diminished for cer- 
tain phase plate configurations [see Fig. [2jc)] and/or by 
employing a large wq (»g ^ V^wrjb)- 

Spiral phase plates were recently used in an elegant ex- 
periment aimed to show the high-dimensional spatial en- 
tanglement of a two-photon state from a down-conversion 
source [14 . The relevant parameters of the experiment 
are: rj = 3.5 for the phase plates, A = 0.8 jiui, pump 
width wo = 780 /im, and thickness of the nonlinear crys- 
tal L = 1 mm |14j . According to the model presented 



in Ref. 



which corresponds to our V 



(SPP) 



it was 



concluded that K. > 3700 ±100. Figure H depicts the 

(SPP) 

predicted distribution for V s (dashed line) together 
with the measured coincidences reported in Ref. [2]. 
According to our approach, a probability distribution 
j-,(spp) g^-^ e( j -^th s = 0.66 (solid line in Fig. [4} shows 
excellent agreement with the experimental results. Equa- 
tion Q then suggests that the corresponding Schmidt 
number should be much smaller than the above K,. How- 
ever, we cannot perform a definite estimate for K, since 
the value of the fiber mode width wq at the location of 
the spiral phase plates, needed to calculate 11 — wq/wg 
appearing in Eq. d7j), was not provided in Ref. |14j . 
The inset in Fig. |4| plots K, for a wide range of ratios 
ji = wq/wq- For instance, using our fitting parameter 
s = 0.66, together with wq — 780 /im and assuming 
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W G > 250 fim, leads to \i < 3 and an estimate for /C < 20, 
at least two orders of magnitude smaller than the quoted 
JC in Ref. [2]. A very large value for JC could only be ex- 
pected if, in the experiment, the fiber mode width icg in 
the plane containing the phase plates was much smaller 
than the pump width wq. 



VI. CONCLUSIONS 

In view of the situation described above, it is rea- 
sonable to conclude that additional experimental results 
and further theoretical analyses -a recent example can 
be found in Ref. |32j - are necessary and of interest in 
order to clarify the problem of quantifying the amount 
of transverse entanglement of photon pairs produced in 
down-conversion sources. If future values for the Schmidt 
number coming from accurate and solid measurements 
confirm the estimate of Ref. |14j and cannot be repro- 
duced within our approach, one should conclude that the 
conventional Gaussian profile constitutes a poor approx- 
imation for the phase-matching function G in Eq. Q. 
However, the ability of our predicted coincidence distri- 
butions to fit the experimental distributions of Ref. [2] 
suggests the correctness of our analysis and that the value 
for the Schmidt number quoted in Ref. [14] could be over- 
estimated. In this case, our approach not only would pro- 
vide a simple an accurate procedure to identify the rele- 
vant spatial modes and extract the degree of transverse 
entanglement; it would go in fact beyond by character- 
izing the action of all local bipartite azimuthal optical 
transformations on a broad family of two-photon states 
using only two detectors. Our results could also be of in- 
terest in other fields such as in identifying the intervening 
spatial modes when preparing well-controlled superposi- 
tions of photon states carrying OAM (to be coherently 
transferred onto Bose-Einstein condensates [33 , 34J or by 
entangling photons with ensembles of cold atoms |35j). in 
high-resolution ghost diffraction experiments with ther- 
mal light [36], or in decoherence processes such as the en- 
tanglement sudden death mediated by the simultaneous 
action of several weak noise sources on bipartite photon 
systems [57] . 
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APPENDIX A 

In this Appendix we outline the derivation of the eigen- 
values and eigenfunctions in the Schmidt decomposi- 
tion pi) when the two-photon amplitude $ is given by 
Eq. (WT. We start with some well-known facts. 

For a general two-photon pure state, \ip} 7 it is possible 
to express \ip) as a bilinear sum of idler and signal basis 
states \Ti :S ) belonging to the Hilbert space of the system 



\ip) = C -rt,r.\n) ® |r s 



(Al) 



where t\ and t s label the set of quantum numbers for 
the idler and signal photons, respectively, whereas the 
coefficients C Ti , Ts describe the probability amplitudes for 
each tensor product of basis states. Our first aim is to 
evaluate the coefficients C TijTa for the down-converted 
state IV') = Jdq i dq s $(q i ,q s )a t (q i )a t (q s )|vac). Let 
u Ti 3 (qi jS ) = (vac|a(qj !S ) \T i s ) denote the basis wave func- 
tions in the transverse momentum representation. Insert- 
ing the closure relations |rj. s )(rj, s | = 1, it can be 
readily seen that the coefficients C TuTs are given by 



/'dq i dq s $(q i ,q s )it*.(q i )u* s (q s ). 



(A2) 



Since we are interested in elucidating the correlation 
properties of the spatial degrees of freedom (OAM and 
the magnitude of the transverse radial wave vectors) , we 
choose as the computational basis the complete set of 
normalised Laguerre-Gaussian modes [2Uj 



w 2 n\ ( wq\ i £ | I w 2 q 2 



_ 27r(|*|+n)! 

/ w 2 q 2 
x exp — I exp 



it<t>- i- (2n + 



exp(i€0)ui j „(g)/v / 27r, 



(A3) 



where q and 4> denote the radial and azimuthal variables 
in momentum space, w is the mode width (at the beam 

waist), and L„ (i) are the associated Laguerre polynomi- 
als. The indices I = 0, ±1, ±2, . . . and n = 0, 1, 2, . . . rep- 
resent the winding (or topological charge) and the num- 
ber of nonaxial radial nodes of the modes. 

Combining Eqs. Q, ( |A2| and {A3} , it follows that 



C££- = ^J dq j; dq se -^l*+^l 2 / 4 e - b2 l*-^l 2 / 4 

x e- u ^e- u «4>°v iuni {qi)vi s ,n a {<ls)- (A4) 

By means of the well-known Anger-Jacobi identity 
e _ a;cos(</>i _ 0s) = ^ = _ oo (-l)" l J m (ar)e iro ^*-^), where 
I m (x) is the modified Bessel function of the first kind, the 
two angular integrals yield the selection rule m — £i — 
—£ s = £. This shows that the idler and signal photons are 
perfectly anticorrelated with respect to their topological 
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charge, which is a manifestation of OAM entanglement. 
Hence, C™'^™ 3 = C""" s Siii S-u s , where 



c?tt = {-iy2w b I I % 

Jo Jo 

ve,m w) v ~t,n s (q s ) It I g J ■ ( 5 ) 



2 _ J.2\ 



The first of the radial integrals in ( A5 ) can be performed 



by resorting to the following formula 

i 

X W +1 LW(x 2 )e-^ 2 I e (xy)dx 

y 2 



L - LJ n 



2\e\+ij\e\+n+i 



4(7-1)7 



valid for all real y, integers n > and £, and complex 7 
(Re( 7 ) > 0). 

The second radial integral can also be done analyti- 
cally. However, a dramatic simplification occurs if the 
widths Wi and w s of the idler and signal radial modes, 
which at this stage have not been specified, are properly 
selected: to, = w s = \f2wob. In this case, one obtains 
a second selection rule for the radial indices. Namely, 
7i{ = n s = n. That is, with such a choice of the widths, 
which is unique, the idler and signal radial modes are 
perfectly correlated. But this shows that it is precisely 
for Wi — w s = \/2w~o~b when one derives the Schmidt de- 
composition. Thus, the Schmidt width, iu<?, corresponds 
to ws — y/2wob. The coefficients reduce then to 



CI 



i-iy 



4wob ( wq — b 



(wo + b) 2 \ Wq + b 



\£\+2n 



(A6) 



Let £ = (wq — b)/(wo+b). It is now clear from Eq. ( A6 ) 
that the Schmidt eigenvalues are X ln = (1 _£2)2£2|*|W 
We therefore conclude that the Schmidt decomposition 
of the two-photon amplitude Q is 



$( qi , qs )=(i-a y, E(- i )^ l+2n ^«(*) 

e=-ac n=0 

x u_4 n (q a ), 



(A7) 



where u*, n (qi) = e I ^w^ ! „(g i )/v27r and u_^„(q s ) = 
e - ll <t>s v _ t n (g s )/y / 27r represent the idler and signal 
Schmidt eigenmodes. These eigenmodes belong to the 
Laguerre-Gaussian basis. It is worth mentioning that 
although the Schmidt width is unique, other decom- 
positions (actually infinitely many) are possible in dif- 
ferent eigenmode bases. For instance, the Hermite- 
Gaussian modes (in Cartesian variables and the same 
ws = \f2w$b) constitute another possible basis. In fact, 
all eigenmode bases connected via the following unitary 
(metaplectic) transformation U(9,<p) = exp(—i6C ■ u v ), 
where = (— sin tp, cos tp, 0) is a unit vector in the 
equatorial plane of the so-called orbital Poincare sphere 



(parameterised by the spherical angles 9 and tp) and C 
is an angular momentum operator, form Schmidt bases 
for ( A7 1 . The components of the angular momentum op- 
erator: C x , C y , and C z are the three SU(2) generators. 
Of these, only C z represents real spatial rotations along 
the propagation of light, and it is thus the only compo- 
nent associated with the orbital angular momentum of 
light that can be measured in experiments |20j . 



APPENDIX B 

In this second Appendix we explicitly show how to cal- 
culate the radial functions TZg(s) given in Eq. (5). For 
the chosen two-photon amplitude Q and the measure- 
ment scenario, these radial functions depend on the spa- 
tial overlap of the Schmidt and the fiber (fundamental 
Gaussian) modes at the planes where the azimuthal phase 
plates are located. We need to evaluate the integrals 

ft< = E n y^Jnr s dndr s v^Zl,n A =o( r i)4^=o,n B =o( r s) 
xv ( W n\' r i) v ^-i S l( r s)- We use the Schmidt eigenvalues 
A £n = (1 - £2)2£2K|+4n found in Appendix A. The idler 

and signal Schmidt modes v^f'^i) an< ^ v -in( r s) be- 
long to the Laguerre-Gaussian mode basis. We stress 
once again the fact that the corresponding fiber rag and 
Schmidt ws widths are not necessarily equal. This im- 
plies that the orthogonality relation for these modes, 



J rdr v 



(r) = S^ejr, 



does not hold in 



general when ws wq- 

We start by recalling a remarkable identity between 

|£| 

Laguerre polynomials L„ and the modified Bessel func- 
tions I 1 [38] 



n\Ll\x)Ll\y) 



(xyz) 



-14/2 



r( 

n=0 y 



exp 



x L 




(x + y)z 



(A8) 



valid for all real x, y, integer £ and complex z [\z\ < 1). 
In detail, the integrals to calculate read 



R< = (i-{ i )f w Es 



w G 



ws 



Kl 



i)3 



2rf\ 2e w s 



irz 



ws Jo 



£\ + n)\ 

r 2 

WG 



+ n)\ 



V2r s 
ws 



1 W 



L W 



2ri\ 2e a s 



w s 



(A9) 



Inserting identity ( A8 1 in ( A9 1 and employing the follow- 
ing result [38] 



re-^h{2 V r)dr = 1 ± 2 f M_ k J- 
2T(l + \£\)Vy/i 2- 2 \ 7 
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valid for Re(7) > and all complex v, with M a< b repre- 
senting the Whittaker functions, we obtain 



Hi = 



r(i + f) 
r(i + \e\) 



dr exp 



8(1 -<e 2 ) 3 - 



(2 - s 2 )£r 2 



)swt 



x M i \i\ 



2s£r 2 



(1-eWs 



(A10) 



where we have introduced the parameter s — 2£/[l+£ 2 + 
(1 — £ 2 )(ws/wg) 2 }- The remaining integral is carried out 
through the use of formula 



< ■ r(i + f) 



oo _ 7 i 
Vt 2,2 



(7-|) (|+7) 



1 



2 2 5 +7 



valid when 7 > 1/2. The function F(a,b; c; d) denotes 
the hypergeometric function. Equation (AlOl is finally 
given by 



Hi = 



i-e 



r 2 (i + f ) 



1 _i_ ( 1- £ 2 ) ( tug. _ wa ] 

Ml 14 

2 ' 2 



r(i 



> r I ! + |*|;» 2 ]« |/| 



(All) 



This expression reduces to Eq. (5l when one leaves aside 
all the ^-independent factors (they do not play any sig- 
nificant role in the normalised coincidence profiles). In 
this way the two-photon detection probabilities ([3| can 
be cast in terms of the impulse response functions t , 

_ e and the radial functions IZe(s) that solely depend 
on the parameter s (when — £g = 0). 



APPENDIX C 

In this appendix we would like to justify the model pro- 
posed here, the double gaussian profile as phase matching 
function. 



1. Classical Energy for Non-linear Media 

In this section we derive explicitly the energy of an 
anisotropic and non-linear medium. This classical de- 
scription enables us to find a reasonable quantized Hamil- 
tonian which allow us to describe the PDC process. In 
anisotropic and non-linear media we must be cautious in 
the description of physical variables. For a deep analysis, 



we begin by assuming that the polarization of the mate- 
rial obeys the usual expansion from non-linear optics 



Pi = eobctfEj + X\ %E 3 E k + xfmEjEkEt 



(A12) 



where (with n > 1) represent the nonlinear suscep- 
tibility tensors of order n+1 responsible for the coupling 
of n + 1 fields. The energy of the electromagnetic field 
can be written as 



W K = d 



E dD , 



(A13) 



where the displacement vector field is defined as usual 
D = e D E + P. In linear media, the energy density of 
the electromagnetic field depends quadratically on the 
electric field cx E 2 , whereas in non- linear (parametric) 
media, polarization exhibits a dependence on higher or- 
der electric field terms (A12). We can expand the energy 



( A13 1 in two terms, one related to the linear contribution 



of the electric field and the other one related to the non- 
linear part. Taking the lowest order nonlinearity ( A12 1 



the energy due to this contribution can be written as 

133 mi, 



Hi 



d 3 r xfjl E l (r,t)E J (r,t)E k (r,t) 



(A14) 



where the integration extends over the volume V of the 
nonlinear medium. Notice that the second order suscep- 
tibility is not the same as the susceptibility in equation 
( A12 1. We will consider that the medium response to 



the electric field is not instantaneous. Hence, the second 
order polarization should be written as 

/oo /*oo 
/ dt 1 dt 2X ( i%(t-ti,t-t2) 
-co J — CO 

x E^v^EkirM) ■ 

(A15) 

We can see, taking the Four ier transformation of the elec- 
tric fields in equation ( |A15 1 , that this specific form of the 
susceptibility implies that different monochromatic waves 
(with a well-defined frequency) of the electric field can 
excite other frequencies. Hence, the energy of the elec- 
tromagnetic field due to the bilinear susceptibility (A14| 
can be expressed as [4Tj . 



Hi 



d A r 



x E i {T,Ly , )E j {T,u-u/)E k {T,u/) . 



(A16) 

This energy gives rise to the quantized Hamiltonian fre- 
quently found in the literature [15] to describe the pro- 
cess of spontaneous parametric down conversion, where 
an incident beam interacts with a non-linear medium and 
splits into two lower-frequency signal and idler photons. 
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2. Parametric Down-Conversion 

If we have a crystal with nonzero x^ 2 \ pumped with a 
laser beam, there is a small, but non negligible, probabil- 
ity that a pump photon will decay into an idler and signal 
photon pair. We take the incident pump light beam to 
be in the form of an intense quasi-monochromatic plane 
wave propagating along the z-direction, 



E„(r,t) 



(A17) 



where e p is the polarization of the pump beam and 
Mnm(ri,z) is a slowly-variant-amplitude (for the mode 
having indices n and m) along z which obeys the paraxial 
equation in anisotropic media [42j . That is, the paraxial 
equation for free propagation cannot be used and must be 
modified. The pump beam is described classically, within 
the undepleted approximation, which means that one as- 
sumes that the input state is made up of a large number 
of photons and only few of them interact with the non- 
linear medium, producing pairs of correlated photons. In 
other words, the pump beam, at the output of the crystal, 
remains almost unaffected. 

A form for the interaction Hamiltonian operator H[{t) 
can be given motivated by the classical electromagnetic 
field energy ( A16 1, and reads (we are always in the inter- 



action picture, though it can also be done in the Heisen- 
berg picture [40]). 



H,(t) 



kj,(Ji k s ,<7 s 



x (4 s , CT J J ( e p) fc Unm(r±,z) al. ^al 



i(q;+q 3 )-rj 



-i[(ki z -\-k sz )z-\-(u! p —uJi—uj s )t] 



B.C. 



(A18) 



where k^ s ) 



,(i,s) 



ki l ' s ^u z . In equation (|A18| two 



quantized electric fields are considered; they will be re- 
sponsible for idler and signal photons, whereas the pump 
beam is associated with the classical field [33] ■ As we 
have mentioned before, we are carrying out the calcula- 
tions under the interaction picture, where we know that 



a state evolves as 

+ GD j {lo dr ll dT ' ^ i{ - t) " i{t ' ) \ | * o>+ -' 

(A19) 

with T denoting temporal ordering. Since the interaction 
is assumed t o be weak enough, it suffices to expand the 
exponential (Al9l up to first order term. The action of 



the Hamiltonian operator on the evolution of the initial 
vacuum state yields the initial state plus a two-photon 
state (contributions due to higher order processes involve 
the generation of four, six,... photons, and are exceed- 
ingly small). Notice that applying the Wick theorem in 
the above expansion, temporal ordering becomes normal 
ordering up to first order. Therefore, 



|*(t))e =T{e-*^ dT ^ (T) } |fi) ~ |0) 



i 1 

hi? 



E E 

k; ,(7^ k s ,(T S 



-i(ujp — 



,)t/a sin [( u p — Ui— Lo s )t/2] 

(iOp ~ LUi — LO s )/2 



ki,<Ti;k s ,CT s ) , 



(A20) 



where we have performed the integration with respect to 
time r. Let L denote the length of the crystal along z, 
whereas S corresponds to the transverse area (the vol- 
ume of the crystal is V = LS) . Equation ( A20 ) is a gen- 
eral expression for any input paraxial beam. Assuming 
that the transverse area of the crystal S is much larger 
than the characteristic width of the input beam, the inte- 
gral with respect to the transverse position becomes the 
transverse Fourier transform. Now taking into account 
that the paraxial wave is propagating in an uniaxially 
anisotropic crystal 02], this integral becomes 



where, for simplicity, we assume that the optical axis 
is contained in the plane xz and the pump beam is lin- 
early polarized along the ir-direction (the pump beam has 
extraordinary polarization). Here, 9 refers to the angle 
between the optical axis and the z-axis, and 



kp — 

c 



eg = cos t 
1 



\{6,ui v ) n\{w p ) n 2 (u p ) h 2 (ujp) n 2 {u p ) n 2 (ui p ) 



Knowledge of the pump spectrum at the input crystal 
face enables, via equation ( A21 1 , to describe the evolu- 



tion of the pump beam through the crystal. In order to 
satisfy both the undepleted approximation and the fact 
that the pump profile is assumed not to vary significantly 
during the nonlinear propagation, the crystal must be 
very thin. The z-dependence in equation ( |A21 1 will play 
a crucial role to correctly describe the phase-matching 
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conditions below. 
Defining 

q_ 1 n e (8,LO p ) / 2 2 n e(^i ti, p) 

2k p nl(u p ) \ y Qx n'l(uj p ) 

we ca n perform the integral with respect to z in equation 
( A20 1, resulting in the so-called phase- matching function: 



h 2 (ujp) 



where the argument 



A(fcp,k i; k s ) = k p ri e (9,u> p ) - 6 - k iz - k sz , (A23) 



jA(fcp,ki,k s )L/2 sin [A(fc p ,kj,k s )L/2] 
A(fc p ,k !; ,k s )/2 



describes the wave vector mismatching. Therefore, equa- 
(A22) tion ( A20 1 can be cast (excluding the vacuum state) as 



x g^q, + q s ) e -*»-H-».W* ^13 | ki ,a i; k s ,a s ). (A24) 

(Up — LUi — U! s )/Z 



3. PDC in type II 

In this section we are going to consider the paramet- 
ric down conversion under type II phase-matching con- 
ditions. This implies that the pump beam is extraor- 
dinarily (linearly) polarized, whereas the idler and idler 
photons are ordinarily (the polarization vector is perpen- 
dicular to the optical axis) and extraordinarily polarized, 
respectively. Our first aim will be to further simplify the 
down-converted state (A24l. If the interaction time is 



sufficiently large, we can replace the temporal-sinc factor 
in equation (A24| by a delta function, 



1 [((jJp - LOj - LO s )t/2] 

(uj p — Ui — uj s )/2 



2tt 5{ujp 



This is a good approximation, and the oscillatory factor 
e -i(wp-wi-w,)t/2 can then discarded by a regulariz- 
ing procedure. In particular, one obtains the frequency 
matching condition 



(A25) 



that is, energy conservation imposes that the annihila- 
tion of the pump photon gives rise to the creation of the 
idler and signal photons. In what follows, and effective 
nonlinear susceptibility Xeff will be used to represent 



(2) 

the relevant components of Xijk contracted with the cor- 
responding components of the pump, idler and signal po- 
larizations (we also absorb in Xef / additional constants) . 
Considering the continuity of the wave vectors of signal 
and idler photons (assuming that the volume of the crys- 
tal is large enough), equation (A24l reduces to 



(2ph) 



1 / j3i / j3i i \ iA(fc„,k, -M,)L/2 

•- / a ki / a k s Xeff(<*> p , m, w s J e ^ *" " s ' ' 

sin[A(fc p ,k„k s )L/2] ra . +Ql 
A(fe PJ k i ,k s )X/2 w "^ +t ^ 



X 5(uip — LOi 



ki,a ,k s ,(j e } 



Notice that the idler and signal photons have fixed po- 
larization, ordinary (<r ) and extraordinary (cr e ) respec- 
tively. 

For convenience, we change the integral in the z com- 
ponent of the wave vector by the frequency. Therefore, 
for the idler integral, the Jacobian due to the change of 
variables is 



d <\idk iz — 



9(q<, k iz 



d <\idui — 



dk l2 



du)i 



d <\idu)i 



u>i n 2 (uji) {J*n (uJi) dn {u>i) 



+ 



ki z C^kiz diOi 

while for the signal integral, we have 



d (\iduji 



2k 2 gz dn c {9,, 



k sz qax sin 26 



d q s dk sz 
( 1 ^ - o 2 

\A 2 (lo 3 ) J ^sx dui 



d(q s , ksz) 
d(q s ,u s ) 



d q s duj s 



dk s 
du), 



i 20 



'nl(0,u> s ) ^ n 2 (w s )^s 



d q s dui s 



d q s dio s 



(A26) 
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Proof: For ordinary waves the dispersion relation in an 
anisotropic medium is 

^nljuji) 2 2 
Solving the quadratic equation for the variable ki z 
h — 

V c z 

Fi- 



Notice that we only choose the positive solution of ki 
nally, taking the derivative with respect to u>i 



dki z 
duJi 



uji n 2 (uji) Lofn (c0i) dn {wi) 



c 2 ki z 



dun 



For extraordinary waves the dispersion relation in an 
anisotropic uniaxial medium is more complex: 

u>s k1 z , . . k az sin 29 

# = nm^j + {qsx cos v + qsv sm V) + 

(q sx cos p + q ay simp) 2 (q ax sin (p - q sy cos ip) 2 . . 

H ~, — s A A ^'J 



n 2 {u a ) 



where {9, tp} are the polar and the azimuthal angles respec- 
tively in spherical coordinates, with 9 referring to the angle 
subtended by the optical axis and the z-axis. For simplicity, 
as we have mentioned before, we assume the optical axis to 
be contained in the plane x-z {tp — 0). We have to take into 
account the common definitions in anisotropic media 









1 


cos 2 6 


sin 2 6» 


n' 2 (6,iu a ) 


n 2 (u! a ) 


n 2 (uj a ) 


1 


1 


1 


n 2 (u) p ) 




n%(u>p) 



In this case, performing the im plicit derivative with respect 
to u) a in the dispersion relation (A27l, 



2^ = 



dk 3 



2k 2 z dn e (9,LO s ) 



n 2 (9,uj a ) du s 



sin 26 dk S z 
Qsx ^57 — r r- k az q sx sin 26 



nf (0, uj s 
d 
duj s 



2 d 

Qsx 



2 
lay 



OuJa 
1 

n 2 (oj s ) 
dn e (uj. 



(A28) 



where 
d 

OuJa 
d 



1 



duj a \n' 2 (9,LU 3 ) J n%{u s ) du a 

1 dn (uj 3 ) 1 dn e (uj s ) 



n 2 (u) s ) 
1 



nl(uja) du>s 



n%(u> s ) duo s 



duia \n2{9,uJs) 



= -2 



cos 2 9 dn e (uja) _j_ sin 2 9 dn (u) s ) 



11% {u)a ) du! s 



nZ(<jj s ) duo s 



Solving the equation ( A28 1 for we find equation 



(A26I 



Substituting equations (A26l into the output state of 
idler and signal photons ( A26 ), we obtain a general ex- 



pression for PDC in type II. So far, we have only assumed 
the undepleted approximation and that the transverse 
area of the crystal is much larger than the characteristic 
width of the pump beam. 



4. PDC in type II, collinear regime 

We now focus on the collinear regime in PDC where 
we only consider those idler and signal photons which 
are produced nearly in the same direction as that of the 
pump beam (in our framework along the z-axis). In this 

particular case, the approximation q \ s L <SC 1 is justified. 



Therefore, equations (|A26|) reduce to 
drq^dki 



for the idler, and 



2 n e (6,uj s ) 2 , 

d q s dk sz = d q s duj s , 

c 

for the signal, where we have also made the reasonable 
assumption that the variation of refractive indices with 
respect to frequency is small. Thus, Jacobians due to 
the change of variables are just refractive indices, which 
depend on frequency. 

In the collinear regime more simplifications can be 
done. Let us explic itly w rite the argument of the phase 
matching function ( A23 1 as 



A(fe p ,ki,k s ) = k p n e (9,Lj p ) 



sin 26> 



nl(9,u v ) nl{9,u a ) 



+ 



k s n e (9,u! s ) — kin (u)i) + 
sin 29 n 2 (9, u> p ) 



+ qix- 



n e (9,uj a ) 
2rii{uj v )k p 

n%(0,u. 



n 2 (uua) 
(q iy + q sy )' 2 + (qix + q S x) 
n e (9, io 



2 n 2 (cj p ) 
2 n 2 e (8,ujp) 



2n 2 (uj s )n 2 (u)a)k 1 



2 

Qax + 



2n 2 (uJa)k i 



2 . 9 
Qsy + 



nl{up) 

2 



+ 



2n (uJi)ki 



(A29) 



It is clearly an strongly anisotropic function. Notice that 
the linear terms in the transverse wave vectors are re- 
sponsible for the walk-off effect (the Poynting and wave 
vectors of the signal and pump photons are not collinear, 
i.e., their energy flows in a different direction from its cor- 
responding phase front). We should mention a previous 
study by Torres and coworkers |43j where the walk-off ef- 
fect was taken into account to describe entangled photon 
states generated in type I PDC. 

The stat e for idler and signal photons within type II 
PDC dA26j) is 



\®}¥ Ph) = / / d 2 <\id0Ji I I d 2 (\aduJa Xeff{Up,LJi,U s ) 

J An Jo J An Jo 

X e iA(fc "' k " ks)L/2 sinc[A(fc p ,k I ,k s )L/2] S„ m ( qi +q s ) 

X S(Wp — LOi — LOa) |ki,<T o ,k s ,0- e ) 

(A30) 

where sinc(a;) = sin(a;)/x. The symbol Af2 represents 
the fact that the transverse idler and signal wave vectors 
arc constrained to be very small (this is the condition for 
the collinear regime). This can be achieved experimen- 
tally by using a suitable pin hole placed in the z axis. 
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Now, due to this restriction, only the maximum of the 
sine function contributes to the integral (the side lobes 
being blocked). Furthermore, the action of the pin hole 
truncating the sine function can actually be modeled by 
a simple exponential function. In addition, we can reg- 
ularize the phase factor e' 4k p ' q ' qa ' by just shifting our 
coordinate system (z = 0) to the middle of the crystal, 
see equation (A23l. We will also focus, for simplicity, on 



the frequency degenerated case, where the idler and sig- 
nal photons have the same frequency {uj p — LUi/2 = u> s /2, 
which can be realized by using appropriate narrow-band 
filters in front of the photodetectors). Thus, the two- 
photon state reads 

|$)(^) = X e//K) J d\ t J d\ s e -A(*p,k 4 ,k.)L/2 

x u nm ((\i +q s ) |k ; , cr ,k s , er e ) . 

(A31) 

The effective susceptibility, due to the energy conserva- 
tion and the frequency degeneration in idler and signal 
photons, now depends only on the frequency of the pump 
beam. 

The argument of the phase-matching function G = 
e -A(fep,ki,k a )L/2 j g s ^y[ too complex. Since we wish to 
maintain our framework at a sufficiently simple level, and 
we will pursue a full analytical study, some further (and 
drastic) approximations are required: we take all refrac- 
tive indices (ordinary and extraordinary) to be identical 
and independent of frequency. Hence, the argument of 
the phase matching function reduces to 



A(/cp, kj, k s ) 



1 



p n eff 



(A32) 



which is, in contrast to equation ( A29 1, an isotropic (with 
cylindrical symmetry) function. The effective refractive 
index n e ff can be interpreted as a suitable fitting pa- 
rameter (in principle, it can be quite different from the 
ordinary and extraordinary refractive indices). 

In spite of the strong approximations performed above, 
it is illustrative to determine their goodness by compar- 
ing the sine and the exponential phase-matching func- 
tions when their arguments are given by equations ( A29 ) 



and (A32), respectively. Figure [5] depicts both profiles. 
It can be seen that for sufficiently small transverse wave 
vectors both profiles are coincident (by suitably choosing 
n e ff)] the deviations arise, typically, when the pin hole is 
already truncating both functions. This qualitatively jus- 
tifies the validity of the approximations made above. Ac- 
tually, this approximation was previously used (without 
any justification) by Law and Eberly [17j to study bipar- 
tite transverse entanglement from type II PDC. Further- 
more, in other previous works, the sine phase matching 
function was even replaced by a constant [HI |3~T1 |4"H I45j . 



In summary, equation ( A31 1 finally becomes (simpli- 
fying the notation and absorbing the constant factors) 



d q, 



j2 

a q a 



*(qi,q s ) |q*,q») , (A33) 

where $ = EG is the two-photon amplitude, with 
E(qi + q s ) = u nm (qi + q s ) denoting the transverse 
(spectrum) profile of the pump beam, and G(q^ — 

q s ) = e ik p"eff ' qi qs ' being the phase-matching func- 
tion. Since we will not focus in later sections on the po- 
larization degree of freedom (recall that idler and signal 
photons are assumed to be linearly polarized), we have 
omitted its explicit dependence in ( A33[ ). 

Equation ( A33 1 will be exploited in Q to carry out 
our discussion about the OAM entanglement shared by 
the downconverted idler and signal photons. 




FIG. 5: Comparison between the sine (solid curve) and the ex- 
ponential (dotted curve) phase-matching functions when their 
arguments are given by equations (A29l and ( |A32 1, respec- 
tively. The shadowed region corresponds to the truncation 
introduced by the pin hole. 
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